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ABSTRACT 

We analyse the Las Campanas Redshift Survey using the integrated conditional 
density (or density of neighbors) in volume-limited subsamples up to unprecedented 
scales (200 Mpc/h) in order to determine without ambiguity the behavior of the density 
field. We find that the survey is well described by a fractal up to 20-30 Mpc/h, but 
flattens toward homogeneity at larger scales. Although the data are still insufficient 
to establish with high significance the expected homogeneous behavior, and therefore 
to rule out a fractal trend to larger scales, a fit with a CDM-like spectrum with high 
normalization well represents the data. 

Subject headings: galaxies: clusters: general - large-scale structure of universe 

Following seminal work of Pietronero and coworkers, (see e.g. Pietronero et al. 1996), the 
possibility of a large-scale fractal distribution of the galaxies has been investigated by various 
authors. In the current literature, however, there are several conflicting estimates of the largest 
scale at which the galaxy distribution can be approximated by a fractal, ranging from a few 
Megaparsecs (e.g. Peebles 1993), to 20 Mpc/h (Davis 1996), to 40 Mpc/h (Cappi et al. 1998), up 
to more than 100 Mpc/h (e.g. Pietronero et al. 1996). A fractal distribution with dimension D is 
characterized by the property that the correlation function 

g(r) = l + Z(r) (1) 

decreases as a power law, ~ r 3 ~ D . Consequently, the average density p c of galaxies at distance r 
from another galaxy , or conditional density, also decreases as ~ r 3 ~ D since, by the definition of 
correlation function, p c = po[l + £(r)], where po is the cosmic average density. 

Naturally, one can infer the scale at which fractality gives way to homogeneity from several 
other observations, like the cosmic backgrounds, although the conclusions are bound to be model 
dependent. The availability in recent years of deep redshift surveys allows finally to study the 
matter distribution directly from its primary tracers, the galaxies. The deepest galaxy redshift 
survey so far published is the Las Campanas Redshift Survey (LCRS), Schectman et al. (1996). 
LCRS contains 23,697 galaxies with an average redshift z = 0.1 , distributed over six 1.5°x80° 
slices. In this paper we determine the behavior of g(r), the volume integral of g(r), and of the 
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fractal dimension D(r), in volume-limited subsamples of LCRS up to r =200 Mpc/h, the largest 
scale so far investigated with such a statistics, by making use of purposedly designed cells. This 
scale is more than four times the scale previously reached by Cappi et al. (1998) using SSRS2. 

Let us begin by discussing why the statistics g(r) is particularly convenient for our purposes. 
By far the most popular two-point estimators that have been used to investigate the clustering of 
the galaxies are the correlation function £(r) and the power spectrum P(k), a Fourier conjugate 
pair. However, the simplest statistics to use to determine the fractal properties of a distribution 
is g(r), or g(r) = V^ 1 Jq g(r')dV, in terms of which a fractal is defined. There is also another 
reason to use g(r) to study the galaxy clustering. The spectrum of a finite survey is actually 
the convolution of the true power spectrum with the survey geometry. As a consequence, both 
the shape and the amplitude of the estimated spectrum are different from the true spectrum. In 
particular, the estimation of the average density from the sample itself forces the spectrum to 
vanish for k — > 0, so that the detection of a turnaround in the spectrum is often suspect (see for 
instance Sylos Labini & Amendola 1996). A similar problem occurs with the correlation function: 
the integral constraint (Peebles 1980) forces the correlation function to become negative at some 
scale, distorting its shape. In fact, using a subscript s to denote quantities estimated in a finite 
sample of size R s , and employing the fact that the density p s in a sample around an observer is a 
conditional density, the following relation holds 

{.w-tw-a*-), (2 ) 

which shows explicitely the integral constraint J £ s dV s = 0. This problem is clearly absent from 
the statistics g(r) because g s and g are simply proportional, g s (r) = g(r)/g(R s ) . It follows that, 
contrary to what happens for £ and P(k), the slope of g s is an unbiased estimator of the slope of 
the clustering trend, i.e. of its fractal dimension. In fact, we can evaluate the fractal dimension as 

D = 3 + p°, (3) 
alogr 

and verify that D s = D. Therefore, to test the claim of fractality it is necessary to check whether 
g(r) can be approximated by a power-law, in which range of scales, and with which slope. The 
approach to homogeneity is then characterized by D — > 3, that is, to a flattening of g(r). 

The statistics g(r) is a differential quantity. When applied to surveys with a relatively small 
density, as are the volume- limited slices of Las Campanas, it tends therefore to be very noisy. In 
this case, one can smear out the noise integrating over cells, obtaining the integrated correlation 
g(r). Here an important problem arises. Consider a distribution of particles described by a 
statistically isotropic correlation function £(r). When we evaluate the expected total number of 
neighbors within a distance r, we are performing the integral 

£to = 4 [ r ti(ry 2 dr>. (4) 
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If the cell is not spherical, the integral becomes £ = / ^(j')W(r', 9, <f>)dV, where the window 
function W(r, 9, (f>) is defined to be constant inside the cell and zero outside and is normalized 
to unity. Of course, only the part of the cell completely contained within the survey has to be 
considered. However, if the cells are not spherically symmetric, the value of £(r) depends on the 
exact form of W(r,9,4>), and in general is different from the definition in (Q). Moreover, W(r, 9, (f>) 
may vary from cell to cell. The obvious solution to this problem is to restrict the analysis to 
spherical cells. However, this limits the scale to the largest sphere contained within the survey 
boundaries which, for most surveys, can be very small: e.g., less than 10 Mpc/h by radius for 
LCRS. As a matter of fact, so far all the works who used the g statistics adopted spherical cells, 
thereby limiting the scales to less than 50 Mpc//i. Since at this scale the fractal behavior is more 
or less within the standard description (e.g., CDM), it is crucial to extend the analysis to larger 
scales. The simplest way to do so is to consider radial cells, that is, cells whose window function 
can be factorized in a radial and an angular function, both normalized to unity 

w(r,e i <f>) = Wr(r)w sl (e,4). (5) 

In this case in fact the function £(r) is the same as in spherical cells, since the angular factor can 
be integrated out. The advantage is that one can design a radial cell such as to maximize the scale 
r, still fitting it within the survey. As elementary as it is, this method has never been used in the 
literature on the galaxy clustering. 

Before going to the data analysis, let us estimate the variance of g. Let us first assume that the 
three-point correlation function can be written as (Peebles 1980) $ijk = Q[£ij£jk + £ij£ik + Cife^ifc] ; 
where Q is independent of the spatial coordinates. Then, the variance of g evaluated in N c 
independent cells containing on average Nq galaxies is (Peebles 1980; Amendola 1998) 



a c 



Nq 1 (1 + £) + o 3 -? + Q@ 1 + 2K 2 ) , (6) 



where a 2 = J W\dV\W2dV2£,i2 • The first term in Eq. (|6|) is the Poisson noise. Inserting the power 
spectrum we have 

a 2 = (27T 2 )- 1 / P(k)W c {k)k 2 dk, (7) 



where W c {k) equals W{k) 2 for spherical cells, but has to be evaluated numerically in the more 
general case that we study here (see Amendola 1998). The last term in the expression (^) can be 
written as (Peebles 1980) Ki = / W\dV\W2dVi^{r\)^(T\2)- In the important case in which £ is a 
power law, it can be shown that a very good approximation is Ki — <r 2 £. For instance, if £ ~ r _1 , 
it turns out that K% = 1.04<r 2 £. In the following, we will always approximate K2 in this way. 

Another problem arises in practice, namely that the cells we use are not independent, both 
because they are partially overlapping and because the clustering scale may be larger than 
the distance from cell to cell. The effect of the correlation is to reduce the number N c at the 
denominator in Eq. (@). For instance, if the cells oversample the volume by a factor of two, it 
means that a cell out of two is redundant, and the effective number of cells can be taken as N c /2. 
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In general, the number of effective independent cells may be approximated as N e = min(A r c , V/V c ), 
where V c is the cell volume, although of course even this is an overestimation of the independent 
cells . Naturally, we could use TV— body simulations to estimate the errors including the cell-to-cell 
correlation, but then we should generate a different simulation for any model we want to compare 
with; moreover, any finite N— body will inevitably cut large scale power, which in the case of 
testing fractals is a particularly severe limitation. Some comparisons with N— body shows that 
Eq. @) underestimates the errors only by 30% at most. 

In the case of an exact fractal, the expected value of g inside a spherical region of radius r 
embedded in a larger box of size Rq is (Coleman & Pietronero 1992) 

g(r) = (r/R ) D - 3 . (8) 

It is not difficult to show that, neglecting a 2 and £ with respect to f 2 , i.e. in the limit of Ro ^> r, 
and neglecting the Poisson noise, the variance in a fractal is 
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where 7 = 3 — D, and ^(7) = 72/[(3 — 7) (4 — 7) (6 — 7)2 7 ]. Notice that the relative error on g is 
independent on the scale, as indeed is found numerically (Amici & Montuori 1998). For instance, 
for Q = 1 and D = 2, as some observations suggest, crj. = N ( T 1 (8/5)g 2 . 

An important consequence of Eq. (^) is that the relative error of the conditional density 
measured in a single cell can be very large for a fractal, more than 100%. Then, the average 
density in a sample of galaxies around us, a conditional density, has such a large variance, in a 
fractal, that it gives in practice no information. A further consequence is that the variance of the 
amplitude of the correlation function, and of related quantities as ro and <7g, makes the use of the 
correlation amplitude, as opposed to its slope, useless in the case of fractals. This problem applies 
also to the case in which the density of a sample n(r) is measured as a function of the distance 
from the observer, without averaging over several cells. 

If we define the scale of homogeneity as the scale at which g flattens so that D > 2.9, then 
we can quantify it in any given CDM model. Clearly, this scale will be larger the higher is the 
normalization erg. It turns out that if eg ~ 1.5, as observed for bright galaxies in SSRS2 (Benoist 
et al 1996) the CDM homogeneity scale can be as large as 50 Mpc/h, and reach ~ 70 Mpc/h for 
clusters (erg « 2): therefore, only at scales larger than 50 Mpc/h the gap between the pure fractal 
model and the standard scenarios begins to be significant. 

The g statistics in spherical cells has been applied to several galaxy surveys (Pietronero, 
Montuori & Sylos Labini 1997, Sylos Labini and Montuori 1997, Sylos Labini et al. 1998, Cappi et 
al. 1998). Here we summarize only the results from the deepest of such surveys, SSRS2 (Da Costa 
et al. 1994) . SSRS2 includes 3600 galaxies in 1.13 sr of the southern sky, down to an apparent 
magnitude of 15.5. The results of the analysis in Cappi et al. (1998) indicate that the conditional 
density decreases as r _1 from 1 to 40 Mpc/h in all the volume limited samples considered, 
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implying fractality with D = 2 on these scales, with no indication of flattening. Similar results are 
presented in Sylos Labini and Montuori (f997) for the APM-Stromlo survey and in Sylos Labini 
et al. (1998) for CfA2 and SSRS2 surveys. These are the largest scales that can be probed via 
this method. In Cappi et al. (1998) the conditional density from the observer alone, i.e. from the 
vertex of the sample, has also been studied, in order to extend the range. The results are that the 
sample is fractal up to 40-50 Mpc//i and tends to flatten above this scale. However, the errors, as 
expected, are quite large; within 2a the deepest samples include all values between D = 2.7 and 
D = 4. The errors quoted in Cappi et al. (1998), moreover, do not include the ensemble variance 
that, for fractals, is very large, as already mentioned. 

We applied our method of the radial cells to the Las Campanas redshift survey, LCRS, 
the deepest redshift survey so far studied. LCRS contains fields which includes galaxies with 
magnitude between 16.0 and 17.3 , and fields with limits 15.0 and 17.7. Every field % has associated 
a sampling factor < fi < 1 which is the fraction of the galaxies randomly chosen out of the total 
number in the field within the magnitude limits; the weights 1/fi must be taken into account 
in the statistics. We considered the only slice with all fields of large magnitude range (slice at 
-12°). We evaluated the conditional density integrated in radial cells with a shape and orientation 
such as not to intersect the survey boundaries. We cut the sample into four VL subsamples, 
denoted as VL147-297, VL 190-330, VL280-410, and VL224-437 (see Fig. 1 and Table I) where the 
numbers give the lower and upper cutoff distance (the two cuts are necessary because LCRS has 
two limiting magnitudes). The use of volume- limited samples avoids the uncertainties connected 
to the radial selection functions. The correlation function in the full magnitude-limited sample 
is evaluated in Tucker et al. (1997). The results for g(r) are shown in Fig. 2. The errors are 
calculated from Eq. (g) using two models, the standard CDM with galaxy normalization a$ g = 1.5 
and r = 0.3 , and the pure fractal with D = 2 and Q = 1. For comparison, we evaluated the errors 
from a simulation of a standard CDM model, and found that our theoretical estimate of the errors 
is approximate to better than 30%. It can be seen that there is an approximate D ~ 2 fractality 
on small scales, up to 20 or 30 Mpc//i, just as in the SSRS2 case, followed by a flattening of the 
slope. The results from SSRS2 of Cappi et al. (1998; sample with magnitude cut at -20) obtained 
with full spheres are reported in Fig. 2 along with VL147-297, the LCRS sample with average 
magnitude closer to SSRS2. As we can see, the scales we reach here are the largest scales ever 
reached for the g(r) statistics. In the case of VL280-410, the trend is decreasing, albeit with a 
change in slope, down to more than 100 Mpc//i, while for VL224-437 (the most sparse sample) a 
very noisy flattening is reached already at 40 Mpc/h. In Fig. 3 we plot the fractal dimension D 
as a function of r, evaluated as least square slope of sets of five contiguous points of g(r). The 
tendence to D = 3 is very clear. It is also seen that the dimension is not really constant in any 
range, although the scatter is too large to infer a clear trend. 

From these results we can conclude that there is a tendency to homogenization around 50-100 
Mpc//i, as expected from a CDM model. However, we remark that we did not detect a clear 
homogeneous behavior, that is g{r) = 1, not even at more than 100 Mpc//i, except for the most 
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sparse sample. This leaves again space for a fractal behavior to larger scales, especially in view of 
the large errors associated to a fractal. The dimension would be however closer to 3 at large scales. 

We compared our results to the CDM-like power spectrum 

P(k) = AkT 2 (T, k)G(n m , n A) a 8g , a 8m , a v ) , (10) 

where T 2 (T,k) is the transfer function of Bond & Efstathiou (1984) , G includes the redshift 
correction of Peacock & Dodds (1994) and the non-linear correction of Peacock & Dodds (1996), 
o v is the line-of-sight velocity dispersion, and the subscripts A, g and m refer to the cosmological 
constant, the galaxies and the total matter, respectively. In Fig. 2 we compare LCRS with a 
tentative model with T = 0.24, a 8g = 1.4, assuming , for the other parameters, the values ag m = 1 
, iltot = 1 ■P'm = 0.4, U\ = 0.6, h = 0.6 , a v = 300 km/sec. Finally, we find the best fit to g by 
varying T and a 8g for scales larger than 30 Mpc//i, in order to avoid the non-linear corrections 
at smaller scales. The results are listed in Table I. The average for all four samples is ag g = 1.5, 

r = 0.3. 

To summarize, in this Letter we have applied the technique of radial cells to LCRS, pushing 
the analysis of fractality to 200 Mpc/h. We have shown that a proper treatment of the errors is 
important in order to compare alternative models, and is particularly crucial for a fractal, in which 
the variance tends to be very large. We have also shown that the use of radial cells allows very 
large scales to be probed. The conclusion is that LCRS is well described by a CDM model with a 
high normalization, up to a 8 ~ 1.6 for the deepest, and thus brightest, sample. The trend can be 
also approximated as a D w 2 fractal up to 20-30 Mpc/h, but shows a clear flattening afterward. 

We thank Ruth Durrer, Michael Joyce, Marco Montuori and Francesco Sylos Labini for very 
interesting discussions throughout the preparation of this work. 
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Fig. 1. — Two of the volume limited subsamples of LCRS that we analyse in this work. 

Fig. 2. — The function g(r) in various VL samples of LCRS with the errors expected in a CDM 
model and in a fractal D = 2 model (dotted lines). The thin line is the tentative CDM fit (see 
text). In the upper right panel we also report for comparison g(r) from a sample of SSRS2, from 
Cappi et al. (1998) 

Fig. 3. — The function D{r) for the four samples, obtained as the least square slope of sets of five 
contiguous values of g(r) 
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